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We analyze the possibility that due to their superfluid properties some compact astrophysical 
objects may contain a significant part of their matter in the form of a Bose-Einstein condensate. 
To study the condensate we use the Gross-Pitaevskii equation, with arbitrary non-linearity. By 
introducing the Madelung representation of the wave function, we formulate the dynamics of the 
system in terms of the continuity equation and of the hydrodynamic Euler equations. The non- 
relativistic and Newtonian Bose-Einstein gravitational condensate can be described as a gas, whose 
density and pressure are related by a barotropic equation of state. In the case of a condensate with 
quartic non-linearity, the equation of state is polytropic with index one. In the framework of the 
Thomas-Fermi approximation the structure of the Newtonian gravitational condensate is described 
by the Lane-Emden equation, which can be exactly solved. The case of the rotating condensate is 
also discussed. General relativistic configurations with quartic non-linearity are studied numerically 
with both non-relativistic and relativistic equations of state, and the maximum mass of the stable 
configuration is determined. Condensates with particle masses of the order of two neutron masses 
(Cooper pair) and scattering length of the order of 10 — 20 fm have maximum masses of the order 
of 2Mq, maximum central density of the order of 0.1 — 0.3 x 10^* g/cm^ and minimum radii in the 
range of 10 — 20 km. In this way we obtain a large class of stable astrophysical objects, whose 
basic astrophysical parameters (mass and radius) sensitively depend on the mass of the condensed 
particle, and on the scattering length. We also propose that the recently observed neutron stars 
with masses in the range of 2 — 2.4Mq are Bose-Einstein Condensate stars. 

PACS numbers: 67.85.Jk, 04.40.Dg, 95.30.Cq, 95.30.Sf 
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I. INTRODUCTION 

At very low temperatures, particles in a dilute Bose 
gas can occupy the same quantum ground state, forming 
a Bose-Einstein (BEC) condensate, which appears as a 
sharp peak over a broader distribution in both coordi- 
nates and momentum space. The possibility to obtain 
quantum degenerate gases by a combination of laser and 
evaporative cooling techniques has opened several new 
lines of research, at the border of atomic, statistical and 
condensed matter physics (for recent reviews see [H, 0I)- 

To say that so many particles are in the same quantum 
state is equivalent in saying that these particles display 
the state coherence. That is, BEC is a particular case of 
coherence phenomena, related to the arising state coher- 
ence. As the gas is cooled, the condensation of a large 
fraction of the particle in a gas occurs via a phase tran- 
sition, taking place when the wavelengths of individual 
particles overlap and behave identically. For the transi- 
tion to take place^particles have to be strongly correlated 
with each other [l|, Q . 

For an ensemble of particles in thermodynamic equi- 
librium at temperature T, the thermal energy of a par- 



ticle is given by fc^T, where ks is Boltzmann's con- 
stant. For a particle of mass m, the thermal wavelength 
is At = ^j2'Kh? /mkBT. Particles become correlated 
with each other when their wavelengths overlap, that is, 
the thermal wavelength is greater than the mean inter- 
particles distance I, Xt > I. The average particle number 
n for particles in a volume V, n = N/V, is related to 
the distance I through the relation nl^ = 1. Hence the 
condition At > I can be rewritten as nA^ > 1, which 
yields the inequality 0] 
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Hence a coherent state may develop if the parti- 
cle density is high enough or the temperature is suffi- 
ciently low. An accurate description of the BEC for 
an ideal gas is based on the Bose-Einstein distribution 
f {p) ~ {cxp [(£p — ^) //cbT] — 1}"^, for particles with 
momentum p, energy £p = p^ jlm and chemical poten- 
tial [i. In the thermodynamic limit A^ — >■ oo, 1^ — >■ oo, 
N jV — ^constant, the fraction of particles condensing to 
the state with p = below the condensation tempera- 
ture Tc is TiQ = 1 ~ {T/Tc)'^^'^ , while ng = above the 
condensation temperature. The condensation tempera- 
ture is Te = 27rn2n2/3/mfcBC^/^, where C = 2.612 [|. 
The dynamical process of Bose-Einstein condensation in 
the canonical ensemble (fixed temperature T) has been 
studied in H. 
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A non-ideal, weakly interacting Bose gas also dis- 
plays Bose-Einstein condensation, though particles inter- 
actions deplete the condensate, so that at zero tempera- 
ture the condensate fraction is smaller than unity, uq < 1. 
A system is called weakly interacting if the characteris- 
tic interaction radius Tint is much smaller than the mean 
inter-particles distance I, rint I- This inequality can 
be rewritten equivalently as nrf^^ ^ 1. If this condition 
holds, the system is called dilute 0- 

Superfluid liquids, like ^He, are far from being dilute. 
Nevertheless, one believes that the phenomenon of su- 
perfluidity is related with BEC. The experimental ob- 
servations and the theoretical calculations estimate the 
condensate fraction for superfluid helium at T = to be 
no ~ 0.10. A strongly correlated pair of fcrmions can be 
treated approximately like a boson. This is why the aris- 
ing superfluidity in '^He can be interpreted as the conden- 
sation of coupled fermions. Similarly, superconductivity 
may be described as the condensation of the Cooper pairs 
that are formed by the electrons or the holes [H . 

An ideal system for the experimental observation of the 
BEC condensation is a dilute atomic Bose gas confined in 
a trap and cooled to very low temperatures. BEC were 
first observed in 1995 in dilute alkali gases such as vapors 
of rubidium and sodium. In these experiments, atoms 
were confined in magnetic traps, evaporatively cooled 
down to a fraction of a microkelvin, left to expand by 
switching off the magnetic trap, and subsequently im- 
aged with optical methods. A sharp peak in the velocity 
distribution was observed below a critical temperature, 
indicating that condensation has occurred, with the al- 
kali atoms condensed in the same ground state. Under 
the typical confining conditions of experimental settings, 
BECs are inhomogeneous, and hence condensates arise 
as a narrow peak not only in the momentum space but 
also in the coordinate space 

If considering only two-body, mean field interactions, 
a dilute Bose-Einstein gas near zero temperature can be 
modelled using a cubic non-linear Schrodinger equation 
with an external potential, which is known as the Gross- 
Pitaevskii equation [H . 

The possibility of the Bose-Einstein condensation has 
also been considered in nuclear and quark matter, in the 
framework of the analysis of the BCS-BEC crossover. 
At ultra-high density, matter is expected to form a de- 
generate Fermi gas of quarks in which the Cooper pairs 
of quarks condensate near the Fermi surface (color su- 
perconductor). If the attractive interaction is strong 
enough, at some critical temperature the fermions may 
condense into the bosonic zero mode, forming a Bose- 
Einstein quark condensate Q . The basic concept of the 
BCS-BEC crossover is as follows: As long as the attrac- 
tive interaction between fermions is weak, the system ex- 
hibits the superfluidity characterized by the energy gap 
in the BCS mechanism. On the other hand, if the attrac- 
tive interaction is strong enough, the fermions first form 
bound molecules (bosons). Then they start to condense 
into the bosonic zero mode at some critical temperature. 



These two situations are smoothly connected without a 
phase transition [lo| . 

One of the most striking features of the crossover is 
that the critical temperature in the BEC region is in- 
dependent of the coupling for the attraction between 
fermions. This is because the increase of the coupling 
only affects the internal structure of the bosons, while 
the critical temperature is determined by the boson's ki- 
netic energy. Thus, the critical temperature reaches a 
ceiling for the large coupling as long as the binding ef- 
fect on the boson mass can be neglected. Even in the 
nuclear matter where the interaction is relatively strong, 
the binding energy of the deuteron is much smaller than 
the nucleon mass. This fact allows us to work within a 
non-relativistic framework to describe such a crossover 
[loj . However, in relativistic systems where the binding 
energy cannot be neglected, there could be two crossovers 
in the relativistic fiuids: one is the ordinary BCS-BEC 
crossover, where the critical temperature in the BEC re- 
gion would not plateau because of the relativistic effect, 
and the second is the crossover from the BEC state to a 
relativistic state, the so-called relativistic BEC (RBEC), 
where the critical temperature increases to the order of 
the Fermi energy [Toj . 

In isospin symmetric nuclear matter, neutron-proton 
{np) pairing undergoes a smooth transition leading from 
an assembly of np Cooper pairs at higher densities to a 
gas of Bose-condcnsed deuterons as the nucleon density 
is reduced to an extremely low value. This transition 
may be relevant to supernova matter or for the crust of 
neutron stars [ll| . A mixture of interacting neutral and 
charged Bose condensates, which is supposed to be re- 
alized in the interior of neutron stars in the form of a 
coexistent neutron superfiuid and protonic superconduc- 
tor, was considered in [l^ . 

The possibility of the existence of some Bose conden- 
sates in neutron stars was considered for a long time (see 
Glendenning [l^ for a detailed discussion). The conden- 
sation of negatively charged mesons in neutron star mat- 
ter is favored because such mesons would replace elec- 
trons with very high Fermi momenta. The in-medium 
properties of the K~ mesons may be such that they 
could condense in neutron matter as well. Bose-Einstein 
condensates of kaons /anti-kaons in compact objects were 
discussed recently [ij, EBl • Pion as well as kaon conden- 
sates would have two important effects on neutron stars. 
Firstly, condensates soften the equation of state above 
the critical density for onset of condensation, which re- 
duces the maximal possible neutron star mass. At the 
same time, however, the central stellar density increases, 
due to the softening. Secondly, meson condensates would 
lead to neutrino luminosities which are considerably en- 
hanced over those of normal neutron star matter. This 
would speed up neutron star cooling considerably p^ . 
Another particle which may form a condensate is the H- 
dibaryon, a doubly strange six quark composite with spin 
and isospin zero, and baryon number two. In neutron 
star matter, which may contain a significant fraction of 
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A hyperons, the A's could combine to form H-dibaryons. 
H-matter condensates may thus exist at the center of 
neutron star s UM . Neutrino superfluidity, as suggested 
by Kapusta [la], may also lead to Bose-Einstein conden- 
sation 

Zero spin bosons, described by real or complex scalar 
fields, are the simplest particles which can be considered 
in the framework of quantum field theory and general 
relativity. Real scalar fields have equilibrium configura- 
tion that were discovered by Scidel and Suen [l^ and 
are called oscillatons. They arc globally regular but arc 
fully time dependent. As for their stability, they seem 
to be quite robust as far as numerical evolution is con- 
cerned [l^. The objects which can be formed by scalar 
fields have been investigated in detail by using mainly 
numerical tools [1^ . Complex scalar fields can form sta- 
ble equilibrium configurations called boson stars [U , 
that are globally regular and whose energy density is time 
independent. The possibility that dark matter is in the 
form of a scalar field [23l - l25l | or a Bosc-Einstcin conden- 
sate [2^ has also been investigated extensively. 

Therefore the physical results presented above show 
that the possibility of the existence of a Bosc-Einstcin 
condensate inside compact astrophysical objects or the 
existence of stars formed entirely from a Bose-Einstein 
condensate cannot be excluded a priori. Such a possibil- 
ity has been in fact suggested recently. Wang [23] used 
the Gross-Pitaevskii equation, together with the associ- 
ated energy functional and the Thomas-Fermi approxi- 
mation, to study a cold star composed of a dilute Bose- 
Einstein condensate. For a static star, the exact solution 
for the density distribution was obtained. A number of 
perturbative solutions for the case of a slowly rotating 
star have also been derived. The effect of a scalar dark 
matter background on the equilibrium of degenerate stars 
was studied by Grifols [2^], with a particular focus on 
white dwarfs, and the changes induced in their masses 
and radii. 

A detailed analytical and numerical analysis of the 
Newtonian Bose-Einstein condensate systems was per- 
formed recently in [2^ and [13], respectively. In [2§] 
an approximate analytical expression of the mass-radius 
relation of a Newtonian self-gravitating Bose-Einstein 
condensate with short-range interactions, described by 
the Gross-Pitaevskii-Poisson system, was obtained. For 
repulsive short-range interactions (positive scattering 
lengths), configurations of arbitrary mass do exist, but 
their radius is always larger than a minimum value. For 
attractive short-range interactions (negative scattering 
lengths), equilibrium configurations only exist below a 
maximum mass. The equation of hydrostatic equilibrium 
describing the balance between the gravitational attrac- 
tion and the pressure due to quantum effects and short- 
range interactions (scattering) was numerically solved in 

It is the purpose of the present paper to develop a 
general and systematic formalism for the study of grav- 
itationally bounded Bosc-Einstcin condensates, in both 



Newtonian and general rclativistic situations. Our ap- 
proach is independent of the nature of the condensate. As 
a starting point we generalize the Gross-Pitaevskii equa- 
tion by allowing an arbitrary form of the non-linearity. 
To obtain a transparent description of the physical prop- 
erties of the BECs we introduce the hydrodynamical rep- 
resentation of the wave function, which allows the formu- 
lation of the dynamics of the condensate in terms of the 
continuity and hydrodynamic Euler equations. Hence the 
Bose-Einstein gravitational condensate can be described 
as a gas whose density and pressure are related by a 
barotropic equation of state. In the case of a conden- 
sate with quartic non-linearity, the equation of state of 
the condensate is given by a polytropic equation of state 
with polytropic index n = 1. In the framework of the 
Thomas-Fermi approximation, with the quantum poten- 
tial neglected, the structure of the gravitational BEG is 
described by the Lane-Emden equation, which can be 
solved analytically. Hence the mass and the radius of the 
condensate can be easily obtained. The case of the ro- 
tating Newtonian condensate is also discussed, by using 
the generalized Lanc-Emden equation. 

By using the equation of state corresponding to the 
Bose-Einstein condensates with quartic non-linearity we 
consider the general relativistic properties of conden- 
sate stars, by numerically integrating the structure equa- 
tions (the mass continuity and the Tolman-Oppcnheimcr- 
Volkoff equation) for a static configuration. In our gen- 
eral relativistic study we consider the cases of conden- 
sates described by both non-relativistic and relativistic 
equations of state, respectively. The maximum mass and 
the corresponding radius are obtained numerically. Bose- 
Einstein condensate stars with particle masses of the or- 
der of two neutron masses (Gooper pair) and scattering 
length of the order of 10 — 20 fm have maximum masses 
of the order of 2Mq, maximum central density of the or- 
der of 0.1 — 0.3 X lO^^g/cm'^ and minimum radii in the 
range of 10 — 20 km. 

The present paper is organized as follows. The Gross- 
Pitaevskii equation is written down in Section [Hi The 
hydrodynamical representation for the study of the grav- 
itationally bounded BECs is introduced in Section IIIII 
The static and slowly rotating Newtonian condensates 
are analyzed in Section IIVI The maximum mass of New- 
tonian condensate stars is discussed in Section |Vl The 
properties of the general relativistic static condensates 
with quartic non-linearity are studied in Section IVII for 
both non-relativistic and relativistic equations of state, 
respectively. The astrophysical implications of our re- 
sults are considered in Section IVlIl We discuss and con- 
clude our results in Section rVIIII 



II. THE GROSS-PITAEVSKII EQUATION FOR 
THE BOSE-EINSTEIN CONDENSATE STARS 

In a quantum system of N interacting condensed 
bosons most of the bosons lie in the same single-particle 
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quantum state. The many-body Hamiltonian describing 
the interacting bosons confined by an external potential 
Vext is given, in the second quantization, by 



H ^ / dr^+ (r) 



2™ 



+ - J dfdf'i'+ {f) (r') V (f - r') ^ (r) ^ (r') , (2) 

where \l/ (f^ and ^+ (r^ are the boson field operators that 
annihilate and create a particle at the position r, respec- 
tively, and V {f — f") is the two-body interatomic poten- 
tial ■ Vrot {r) is the potential associated to the rotation 
of the condensate, and is given by 



Vrot (r) = frot (t) 



(3) 



where uj is the angular velocity of the condensate and 
frot (t) a function which takes into account the possible 
time variation of the rotation potential. For a system 
consisting of a large number of particles, the calculation 
of the ground state of the system with the direct use of 
Eq. ^ is impracticable, due to the high computational 
cost. 

Therefore the use of some approximate methods can 
lead to a significant simplification of the formalism. One 
such approach is the mean field description of the con- 
densate, which is based on the idea of separating out 
the condensate contribution to the bosonic field opera- 
tor. For a uniform gas in a volume V, BEG occurs in the 
single particle state ^'o = \/ N/V, having zero momen- 
tum. The field operator can then be decomposed in the 
form {r) ~ y/ N/V + ^' (f). By treating the operator 
^' (r) as a small perturbation, one can develop the first 
order theory for the excitations of the interacting Bose 
gases [H, Si- 

In the general case of a non-uniform and time- 
dependent configuration, the field operator in the Heisen- 
berg representation is given by 



^ (r, t) = V (r, t) + (r, t) 



(4) 



where ip{f,t), called the condensate wave function, is 
the expectation value of the field operator, -0 (r, t) = 
(^(f, i)). It is a classical field and its absolute value 
fixes the number density of the condensate through 
n{f,t) = I"!/) (r, The normalization condition is 

N = J n (rjt) dr, where N is the total number of par- 
ticles in the star. 

The equation of motion for the condensate wave func- 
tion is given by the Heisenberg equation corresponding 
to the many-body Hamiltonian given by Eq. ([2|), 



y^^Vrot {r) + Vext + 

Zra 

dr'i!+ {r',t)V{f' -r)^{r',t) 



i>{r,t). (5) 



Replacing ^ (f, t) by the condensate wave function ip 
gives the zeroth-order approximation to the Heisenberg 
equation. In the integral containing the particle-particle 
interaction V {r' — r) this replacement is in general a 
poor approximation for short distances. However, in a 
dilute and cold gas, only binary collisions at low energy 
are relevant and these collisions are characterized by a 
single parameter, the s-wave scattering length, indepen- 
dently of the details of the two-body potential. There- 
fore, one can replace V {•f — r) hy a.'o. effective interaction 
V {f — f) = XS {r' — f), where the coupling constant A is 
related to the scattering length a through A = Anh^a/m. 
Hence, we assume that in a medium composed of scalar 
particles with non-zero mass, the range of Van der Waals- 
type scalar mediated interactions among nucleons be- 
comes infinite, when the medium makes a transition to a 
Bose-Einstein condensed phase. 

With the use of the effective potential the integral in 
the bracket of Eq. (O gives A I?/" (f, and the result- 
ing equation is the Schrodinger equation with a quartic 
nonlinear term [HQ- However, in order to obtain a more 
general description of the Bose-Einstein condensate stars, 
we shall assume an arbitrary non- linear term g' {\il){f,t)\ ) 
[3l| . As was pointed out in [3l|, the Gross-Pitaevskii 
approximation is a long-wavelength theory widely used 
to describe a variety of properties of dilute Bose con- 
densates, but for short-ranged repulsive interactions this 
theory fails in low dimensions, and some essential modi- 
fications of the theory are necessary. 

From a physical point of view these modifications can 
be understood as follows [3l|. The inter-particle inter- 
action can be written as V [f] = uq5^ (r) , where uq is 
the amplitude of the inter-particle repulsion and (r) 
denotes any well localized d-dimensional function that 
transforms into the mathematical Dirac delta distribu- 
tion when the range of interactions a — > 0. Assume 
that the inter-particle interaction is so strong that each 
particle is localized within ge formed by its neigh- 
bors. In the dilute limit nl'^ ^ 1 the size of this cage 
can be estimated as i? n^-^^''- and the ground state 
energy per particle follows from the uncertainty prin- 
ciple as lT?/mR^ ^ f?r?l'^ jm. The ground state en- 
ergy which would go into the energy functional is given 
by W'n^'^^'^^^'^ jm. The strong interaction assumption is 
valid if the interaction energy per particle uq/ R^^ is much 
bigger than the ground state energy per particle, i.e. 
uo/W^ » /inR^. The condition for the strong coupling 
limit can be written as h'^n^'^~'^^^'^ /muQ <^ 1. As space 
dimensionality decreases, it becomes increasingly harder 
for the repulsive particles to avoid collisions. Thus, in 
general, the quartic non-linearity in the energy functional 
should be replaced by . The non-linearity may 

also be logarithmic in [3l|. For d = 1 we have a jV'l^ 
interaction. 

Therefore the generalized Gross-Pitaevskii equation 
describing a gravitationally trapped Bose-Einstein con- 
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densate is given by 



Zm 



(6) 



where we denoted g' = dg/dn. 

As for Vext i'f^, we assume that it is the gravitational 
potential, Vext = and it satisfies the Poisson equa- 
tion 



where 



p = mn = (?% t)\ , 



(7) 



(8) 



is the mass density inside the Bose-Einstein condensate 
star. 



III. THE HYDRODYNAMICAL 
REPRESENTATION OF THE BOSE-EINSTEIN 
GRAVITATIONAL CONDENSATE 

The physical properties of a Bosc-Einstcin condensate 
described by the generalized Gross-Pitaevskii equation 
given by Eq. © can be understood much easily by using 
the so-called Madelung representation of the wave func- 
tion [l], Q , which consists in writing -0 in the form 



■41 (f, t) = a/ n (f , t) 



exp 



(9) 



where the function S (r, t) has the dimension of an action. 
By substituting the above expression of the wave function 
into Eq. ([6]) it decouples into a system of two differential 
equations for the real functions n and v, given by 



— + V • (nt;) = 0, 



(10) 



dv 



+ VQ + Kot + K.t + g' =0, (11) 



(12) 



dt \ 2 

where we have introduced the quantum potential 

2m \fn 

and the velocity of the quantum fluid 
- 



(13) 



respectively. From its definition it follows that the veloc- 
ity field is irrotational, satisfying the condition V x i/ 0. 
The quantum potential Vg has the property [l| 



(14) 



nViVg — \I j ( —-^ — 7T.ViVj Inn j = Vjcr, 



where a^^ = — (?i^n/4m)ViVj Inn is the quantum stress 
tensor, which has the dimension of a pressure and is an in- 
trinsically anisotropic quantum contribution to the equa- 
tions of motion. 

By taking into account that the flow is irrotational, 
the equations of motion of the gravitational ideal Bose- 
Einstein condensate take the form of the equation of con- 
tinuity and of the hydrodynamic Euler equation, 



dp 

at 



^ + V • ipv) = 0, 



(15) 



dv 
di 



+ {v ■ V) V 



Vrot 

m 



ext 

m 



where we have denoted 



my \m/ m \m 



(16) 



(17) 



Therefore the Bose-Einstein gravitational condensate can 
be described as a gas whose density and pressure are 
related by a barotropic equation of state [l| . The explicit 
form of this equation depends on the form of the non- 
linearity term g. 

For a static ideal condensate, v = 0. In this case, from 
Eq. pT|) . we obtain 

Vq + Vrot + Vext + 5' = constant (18) 

Applying the operator to both sides of Eq. (|18p gives 

{Vq + Vrot + 9') + V^Vext = 0. (19) 

In the case of a condensate with a non-linearity of the 
form g{n) = kom^n^ /2, where is a constant, and in 
the presence of a confining gravitational field Vext = m^, 
it follows that the generalized potential Vgen = — Vq — 
Vrot — komp satisfies the Poisson equation. 



1 



gen 



inGp. 



(20) 



If the quantum potential can be neglected, then from 



Eq. by using the relations V-^K-i 



and 



V^g' = komW^p, it follows that the mass density of the 
condensate is described by a Helmholtz type equation, 
given by 



AttG 
ko ' 



UJ 

ko 



0. 



(21) 



IV. STATIC AND SLOWLY ROTATING 
NEWTONIAN BOSE-EINSTEIN CONDENSATE 
STARS 

When the number of particles in the gravitation- 
ally bounded Bose-Einstein condensate becomes large 
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enough, the quantum pressure term makes a significant 
contribution only near the boundary of the condensate. 
Hence it is much smaller than the non-linear interaction 
term. Thus the quantum stress term in the equation 
of motion of the condensate can be neglected. This is 
the Thomas- Fermi approximation, which has been exten- 
sively used for the study of the Bosc-Einstcin condensates 
[l[ . As the number of particles in the condensate becomes 
infinite, the Thomas-Fermi approximation becomes exact 
j27l [29I [soj . This approximation also corresponds to the 
classical limit of the theory (it corresponds to neglecting 
all terms with powers of H) or as the regime of strong 
repulsive interactions among particles. From a mathe- 
matical point of view the Thomas-Fermi approximation 
corresponds to neglecting all terms containing Vp and 
VS* in the equation of motion. 



A. Static Bose-Einstein Condensate stars 

In the case of a static Bose-Einstein condensate, all 
physical quantities are independent of time. Moreover, 
in the first approximation we also neglect the rotation 
of the star, taking Vrot = 0. Therefore the equations 
describing the static Bose-Einstein condensate in a grav- 
itational field with potential Vext = "m-^ take the form 



VP 



(22) 



(23) 



These equations must be integrated together with the 
equation of state P = P (p), which follows from Eq. ([T7| , 
and some appropriately chosen boundary conditions. By 
assuming that the non-linearity in the Gross-Pitacvskii 
equation is of the form 



'in) 



(24) 



where a and 7 are positive constants, it follows that the 
equation of state of the gravitational Bose-Einstein con- 
densate is the polytropic equation of state. 



P{p) =a(7-l)n'^ = Kp'', 



(25) 



Hence the mass and the radius of the condensate are 
given by 



R 



in + l)K' 
AttG 



1/2 



(l-n)/2n 



(27) 



and 



M = 47r 



{n + l)K 
AttG 



-1 3/2 



^(3-„)/2„^2|^,(^^)|^ (28) 



respectively, where defines the zero- pressure and zero- 
density surface of the condensate: 9 {£,1) = [s^ . 

In the standard approach to the Bose-Einstein conden- 
sates, the non-linearity term g is given by 



Uq 



(29) 



where uq = infi^a/m 0]. The corresponding equation of 
state of the condensate is 



with 



27rft2 



Pip)^Kp^ 



0.1856 X 10^ 



1 fm / V 2m„ 



(30) 



(31) 



where m„ = 1.6749 x 10"^'' g is the mass of the neutron^. 

Therefore, the equation of state of the Bose-Einstein 
condensate is a polytrope with index n — 1. In this 
case the solution of the Lane-Emden equation can be 
obtained in an analytical form, and the solution satisfying 
the boundary condition (0) = 1 is 



0(0 = 



sin^ 



(32) 



The radius of the star is defined by the condition 
^(Ci) = 0, giving ^1 = TT. Therefore the radius R of 
the Bose-Einstein condensate is given by 



i? = TT 



Gm? 



6.61 



1 fm 



1/2 



2nir 



-3/2 



km. (33) 



where we denoted iiT = a (7 — 1) /m"' . 

By representing 7 in the form 7 = 1 -I- 1/n, where n is 
the polytropic index, it follows that the structure of the 
static Bose-Einstein condensate star is described by the 
Lane-Emden equation. 



For m = 2m„ and a = 1 fm the radius of the condensate 
is i? w 7 km. The radius of the gravitationally bounded 
Bose-Einstein condensate is independent on the central 
density and on the mass of the star, and depends only on 
the physical characteristics of the condensate. 



1 d 



(26) 



where is a dimensionless variable defined via p = PcO"", 
^ is a dimensionless coordinate introduced via the trans- 
formation r = [{n + 1)K pl^"'~^ /AirGj^^'^^ and pc is the 
central density of the condensate [11] . 



^ In adopting this scaling, wc have in mind the possibility that neu- 
trons in the core of neutron stars form the equivalent of Cooper 
pairs and behave as bosons of mass 2m„. This means that we 
treat the core of neutron stars as a superfluid (see Section IVIIII 
for additional comments). However, our study may be valid in 
other circumstances so that we leave the mass m unspecified. 
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The mass of the star is obtained as 

3/2 



M 



(34) 



where we have used \9' (^i)| = I/tt. As a function of the 
central density and of the coherent scattering length a 
we obtain for the mass of the Bose-Einstein condensate 
star with quartic non-linearity the expression 



M = 1.84 



10 



16 



g cm 



1 fm 



3/2 



2to„ 



-9/2 



Mo 



= 5 X 10^ 



^/cm" 



(35) 
, the 



For m = 2m„, a = 1 fm and 
mass of the condensate is M w 0.92 A/q. However, this 
mass may be larger than the maximum mass allowed by 
general relativity (see below), and for a correct determi- 
nation of the maximum mass of BEC stars we cannot 
ignore the effects induced by the space-time curvature. 

The mass of the static condensate can be expressed in 
terms of the radius and central density by 



TT 



(36) 



which shows that the mean density of the star p = 
3Af /AttR^ can be obtained from the central density of 
the condensate by the relation p = Spc/ir^. 

With respect to a scaling of the parameters m, a and 
Pc of the form m — aim, a — ^ Q!2a, Pc — ^ «3Pc, the 
radius and the mass of the condensate have the following 
scaling properties: 



3/2 1/2 



R, M ->■ 



9/2 3/2 



asM. (37) 



B. Slowly rotating Bose-Einstein Condensate stars 

The case of slowly rotating Bose-Einstein condensates 
can also be straightforwardly analyzed, by taking into 
account the fact that the condensate obeys a polytropic 
equation of state. The study of the slowly rotating poly- 
tropcs was performed in detail in 



The Lane-Emdcn equation for a rotating Bose-Einstein 
condensate is 



1 d 



1 d 



(38) 

where p = cos 6* and = uP' jl-nGpc- The volume Vui and 
the mass M^j of the condensate in slow rotation are given 
in the first order in f2 by 



3^0 (Ci) 
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and 



eii^'ieoi 
gi/3~V'(, (a) 



(39) 



(40) 



respectively. In these expressions, Mq is given by Eq. 
(Pg)) and Vb = (4/3)7ri?3 where i?o is given by Eq. 
The values of the function t/iq are tabulated in [s^ . Equa- 
tions ([M)) and (Uni) represent the mass and volume rela- 
tions for two stars with equal central density, one rotating 
with an angular velocity uj and the other non-rotating. 

In the case of Bosc-Einstcin condensates with quar- 
tic non-linearity, corresponding to a polytropic index 
ri = 1, the Lane-Emden equation can be integrated ex- 
actly (yielding di^) = sin^/^ and ■!/'o(0 = 1 - sin^/^), 
giving for the volume Kj and mass M^j of the rotating 
condensate the following simple relations 



M,, 



v^^Vo{i + 2,n). 



Mo 







- 1^ 17 


1 + 


(f 









(41) 



(42) 



V. MAXIMUM MASS OF THE STATIC 
NEWTONIAN BOSE-EINSTEIN CONDENSATE 
STARS: QUALITATIVE TREATMENT 

The numerical values of the basic parameters (mass 
and radius) of the condensed object sensitively depend 
on the mass m of the particle, on the scattering length 
a, and on the central density pc- R = R{a,m), M = 
M (a, m, Pc)- Of course, in general, the values of the mass 
and radius of the gravitational condensate depend on the 
adopted model for the non-linearity. 

The scattering length a is defined as the zero-energy 
limit of the scattering amplitude /, and it can be related 
to the particle scattering cross section a by the relation 
cr = 47ra^ [l|. On the other hand, the notion that par- 
ticles like, for example, the quarks, retain their usual 
properties and interactions at the very hig h densities in 
the neutron stars may not be viable [l3l Is^ . In our 
calculations we use a "hard" core approximation of the 
potential. Therefore we accept that at high densities the 
"hard" core potential is in the QCD range of 1 fm and the 
allowed values of the scattering length a may generally 
be in the interval 0.5 fm < a < 1 — 2 fm, corresponding 
to a scattering cross section of around 1 mb. 

The transition temperature to a Bose Einstein Con- 
densate of dense matter can be written as 



P' 



,2/3 



P 



= 1.650 X 10' 

2/3 



1016 g/cm3 



2m„ 



-5/3 



K. (43) 



Neutron stars are born with interior temperatures of the 
order of 2 — 5 x 10^^ K, but they rapidly cool down via 
neutrino emission to temperatures of less than 10^° K 
within minutes. Also strange matter, pion condensates, 
A hypcrons, 5 isobars, or free quark matter might form 
under the initial thermal conditions prevailing in the very 
young neutron star. Hence a condensation process can 
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take place in the very early stages of stellar evolution. 
If the core is composed of only "ordinary" matter (neu- 
trons, protons, and electrons), then when the tempera- 
ture drops below about 10^ K all particles are degener- 
ate. We expect that after a hundred years or so the core 
will become superfluid [l^l, and this may also favor the 
possibility of a Bose-Einstein Condensation through the 
BCS-BEC crossover. 

Restriction on the maximum central density and max- 
imum mass of the condensate quartic non-linearity can 
be obtained from the study of the speed of sound, de- 
fined as = dP/ dp. With the use of Eq. ([30|) we obtain 

= 2Kp. The causality condition implies Cg < c, where 
c is the speed of light. 

By introducing the dimensionless parameter k, defined 

as 



1 fm 



L/2 / 



-3/2 



(44) 



With a = 1 fm and m = 2m„ we obtain for the maximum 
mass of the Bosc-Einstcin condensate star the restriction 
M < 2.23Mq. 

For the n — \ polytropc the radius of the star is inde- 
pendent on the central density. Generally, one may con- 
sider a as a free parameter, which must be constrained 
by the physics of the nuclear interactions taking place 
in the system. However, due to the possible dependence 
of the free scattering length a on the mass density, in 
the case of Bose-Einstein condensates there may be (in- 
direct) dependence of the radius on the central density 
of the star. 

Finally, we would like to point out that the estimates 
on the maximum mass obtained in the present Section 
are qualitative with respect to the numerical factors, and 
more precise values of the maximum mass of the BEC 
stars will be obtained in the next Section by using a full 
general relativistic approach. 



the causality condition gives the following upper bound 
for the central density of the condensate: 



Pc < 



9 

m c 



2.42 X 10'" K"^ g/cm^ 



(45) 



With the use of Eqs. (|45l) and ()34p we obtain the following 
restriction on the maximum mass of the Bose-Einstein 
condensate with quartic non-linearity: 



M < TT- 



(Gm) 



3/2 



4.46 K Mr. 



(46) 



For m = 2m„ we obtain the condition pc < 
[2.42/a(fm)] x 10'" g/cm^. By taking into account that 
for this range of high densities a physically reasonable 
value for the scattering length is a w 1 fm, we ob- 
tain the restriction on the maximum mass of the Bose- 
Einstein condensate star from the causality condition as 
M < 4.46M0. 

A stronger bound on the central density can be derived 
from the condition that the radius of the star R must be 
greater than the Schwarzschild radius Rs = 2GAf/c^, 
R > Rs. For a Bose-Einstein condensate star Rs can 
be expressed as a function of the central density and of 
the radius as Rs = SGR'^Pc/ttc'^ . Then the condition of 
stability against gravitational collapse gives 



< 



87raft2 



1.21 X 10 



16 , 



?/cm^ (47) 



a relation which for the condensate star with m = 2m„ 
and a = 1 fm leads to the constraint pc < 1.21 x 10'" 
g/cm^. The constraint on the maximum mass for the 
stellar type Bose-Einstein condensate can be formulated 



M < 



2(Gm)3/2 



2.23 kM, 



0- 



(48) 



VI. GENERAL RELATIVISTIC 
BOSE-EINSTEIN CONDENSATE STARS 

In the previous Sections we have considered the gravi- 
tationally bounded Bose-Einstein condensate stars in the 
framework of Newtonian gravity. General relativistic ef- 
fects may change the physical properties of compact ob- 
jects in both a qualitative and quantitative way. For ex- 
ample, general relativity imposes a strict limit on the 
maximum mass of a stable compact astrophysical ob- 
ject, a feature which is missing for classical Newtonian 
stars. Therefore the study of the general relativistic 
Bose-Einstein condensates offers a better understanding 
of their physical properties. In the present Section, we 
study the properties of static general relativistic Bose- 
Einstein condensate stars. 



A. Static general relativistic BEC stars 

For a static spherically symmetric star, the interior line 
element is given by 



-{r)dt^ + e'^('')dr2 + {de"^ + sin^ Odcj)^) . (49) 



The structure equations describing a general relativis- 
tic compact star are the mass continuity equation and the 
Tolman-Oppenheimer-VolkofF (TOV) equation of stan- 
dard general relativity, and they are given by [l3| 



dM 
dr 



Anpr , 



dP{r) G{p + P/c^) {'i-KPr^/c'^ + M) 



dr 



r2 [1 -2GA'/(r)/c2r] 



(50) 



(51) 



These equations extend the classical condition of hydro- 
static equilibrium for a self-gravitating gas to the context 
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of general relativity. We have written the energy density 
as e = pc^. The system of equations ([5n)) - ([5T|) must be 
closed by choosing the equation of state for the thermo- 
dynamic pressure of the matter inside the star, 



P = P{p) 



(52) 



At the center of the star, the mass must satisfy the 
boundary condition 



Af(0) 0. 



(53) 



For the thermodynamic pressure P we assume that it 
vanishes on the surface, P{R) = 0. 

The exterior of the Bosc-Einstein condensate star is 
characterized by the Schwarzschild metric , describing the 
vacuum outside the star, and given by [l3| 



1 



2GM 



r> R. 



(54) 



The interior solution must match with the exterior so- 
lution on the vacuum boundary of the star. 



B. Mciximum mass of non-relativistic BECs with 
short-range interaction: n — 1 polytropes 

We assume that, in general relativity, the BEC can still 
be described by the non-relativistic equation of state 



P = Kp^ 



with 



K 



(55) 



corresponding to a polytropic equation of state, with 
polytropic index n = 1. The theory of polytropic 
fluid spheres in general relativity has been developed by 
Tooper [35| . and we shall use his formalism and nota- 
tions. Therefore, we set 



PcO, p = Kp^ 



A' 



M{r) 



4:TTPc 
A3 



Kp, 



27rG 



1/2 



(56) 



(57) 



where pc is the central density. In terms of these vari- 
ables, the TOV equation and the mass continuity equa- 
tion become 



(58) 



^'(O) = 0. On the other hand, the density vanishes at 
the first zero of 6': 



0(6) = 0. 



(60) 



This determines the boundary of the sphere. In the non- 
relativistic limit CT — 7> 0, the system (|58 |) -(|59 |) reduces to 
the Lanc-Emden equation (|26|) with n ~ 1. 

From the foregoing relations, we find that the radius, 
the mass and the central density of the configuration are 
given by 



R = Cii?*, M = 2crw(6)A'/*, Pc = o-p*, 



(61) 



where the scaling parameters i?*, il/* and p^ can be ex- 
pressed in terms of the fundamental constants and the 
parameter k as 



Gm? 



1/2 



= 2.106 k km. 



(Gm) 



3/2 



= 1.420 kM, 



= 4.846 X 10^6 g/cm^ 



(62) 



(63) 



(64) 



We note that the expression of the scaled radius i?* is 
the same as in the Newtonian regime (in particular it is 
independent on c) , while the scaling of the mass and of 
the density are due to relativistic effects. By varying a 
from to -|-oo, we obtain the series of equilibria in the 
form A'/(pc), R{pc) and M{R). 

The velocity of sound is = P'{p) = IK p. The con- 
dition that the velocity of sound at the center of the con- 
figuration (where it achieves its largest value) is smaller 
than the velocity of light can be expressed as 2Kpc < c^, 
or equivalently a& a < cis with 



(65) 



The values of ^ and u(6) at this point have been tabu- 
lated by Tooper (and confirmed by our numerical study): 



6 ^ 1.801, w(6) = 0.4981. 



(66) 



The corresponding values of radius, mass and central 
density are 

/ fc2 \ 1/2 

i?s = 1.801 f ^— J j = 3.790 k km, (67) 



(59) 



For a given value of the relativity parameter cr, they 
have to be solved with the initial condition 9{Q) = 1 
and v{0) = 0. Since w ^ as ^ — >■ O, it is clear that 



M, = 0.498 



0.707 kM, 



{Pc)s 



(68) 



2.423 X lO^V"^ g/cm^ (69) 
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However, it is not granted that the criterion a > is 
equivalent to the condition of dynamical instability. The 
principle of causality is a necessary, but not a sufficient, 
condition of stability^. The condition of dynamical insta- 
bility corresponds to the turning point of mass dM = 
and there is no reason why this should be equivalent to 
Cs = c. In fact, our numerical study demonstrates that 
this is not the case. We find that the maximum mass docs 
not exactly correspond to the point where the velocity of 
sound becomes equal to the velocity of light. In the series 
of equilibria (parameterized by the central density ct), the 
instability occurs sooner than predicted by the criterion 
(|65|). We find indeed that instability (corresponding to 
the mass peak) occurs for a > cfc with 

CTc = 0.42. (70) 

The values of and u(^i) at this point are 

^1 = 1.888, v{^i) = 0.5954. (71) 

The corresponding values of radius, mass and central 
density are 

i?mi„ = 1.888 ( ^-^ ) = 3.974 k km, (72) 



0.5 


- 


M 
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Relativistic 
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0.2 








0.1 

n 




Newtonian 





0.001 0.01 0.1 1 10 100 1000 

FIG. 1: Dimensionless mass-central density relation of a rela- 
tivistic BEG with short-range interactions modeled by a n = 1 
polytrope. There exists a maximum mass Mmax/Af, = 0.5001 
(black bullet) at which the system becomes dynamically un- 
stable. This mass does not coincide with the point (white 
bullet) at which the velocity of sound becomes equal to the 
velocity of light (see Fig. In particular, gravitational insta- 
bility occurs slightly sooner than what is predicted from the 
criterion based on the velocity of sound (principle of causal- 

ity). 



M„ 



0.5001 



= 0.710 kMq, (73) 



(Pc)r 



0.42 



27ran2 



= 2.035 X lO''' g/cTx^\ (74) 



We also note that the radius of a BEC star is necessarily 
smaller than 



R-n 



Gm3 



= 6.61 K km. 



(75) 



corresponding to the Newtonian limit {a — 0). There- 
fore, its value 3.974k < i?(km) < 6.61k is very much 
constrained. 

The dimensionless curves giving the mass-central den- 
sity, radius-central density, mass-radius relations and 
some density profiles are plotted in Figs. [T][5l In Fig. [5] 
we present the mass-radius relation for a = 1 fm and 
different values of m. 



C. Maximum mass of relativistic BECs with 
short-range interaction 

The previous treatment is approximate because we use 
the equation of state (f55|) obtained in the non-relativistic 



^ Similarly, in statistical mechanics, the condition that the specific 
heat is positive in the canonical ensemble is a necessary, but not 
a sufficient, condition of canonical stability [ssll . 
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FIG. 2: Zoom of Fig. [T]showing that the maximum mass does 
not coincide with the point at which the velocity of sound be- 
comes equal to the velocity of light. The system is stable (S) 
before the turning point of mass (cr < Gc) and unstable (U) af- 
ter the turning point of mass (a > CTc)- The velocity of sound 
is smaller than the velocity of light before the white bullet 
(cr < (Ts) and larger after that point (a > Gs). Therefore, 
there exists a small region (tJc < a < Gs) where the system 
is unstable although the velocity of sound is smaller than the 
velocity of light. 



regime (i.e. from the Gross-Pitaevskii equation) but solve 
the TOV equation expressing the condition of hydrostatic 
equilibrium in general relativity. A fully relativistic ap- 
proach based on the Klein-Gordon-Einstein system has 
been developed by Colpi et al. [l^ . In order to make the 
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FIG. 3: Dimensionless radius-central density relation of a 
relativistic BEG with short-range interactions modeled by a 
n = 1 poly trope. 
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FIG. 5: Dimensionless density profiles corresponding to cr = 
(Newtonian), a = ctc = 0.42 (maximum mass) and cr = (Js = 
1/2 (where Cs = c). 



a = 1 fm 




8 10 
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FIG. 4: Dimensionless mass-radius relation of a relativistic 
BEG with short-range interactions modeled by a n = 1 poly- 
trope. There exists a maximum mass Mmax/M, = 0.5001 
and a minimum radius Rmi-n/ R* ~ 1.888 corresponding to a 
maximum central density (pc)max = 0.42p,. There also ex- 
ists a maximum radius i?max/i?* ~ tt corresponding to the 
Newtonian limit cr — !> 0. 



FIG. 6: Mass-radius dependence for general relativistic Bose- 
Einstein condensates with quartic non-linearity for a = 1 fm 
and different values of the mass m. From top to bottom: m = 
m„, m = 1.25m„, m — 1.5mn, m = 1.75m„ and m = 2m„. 
For all configurations pc > Pn, where p„ — 2.026 x 10^* g/cm'^ 
is the nuclear density, and the causality condition < c is 
satisfied. 



correspondence between BECs with short-range inter- 
actions described by the Gross-Pitaevskii equation and 
scalar fields with a ■jA|0|"' interaction described by the 
Klein-Gordon equation, we set [1^ 



a 



(76) 



where Ac = h/mc is the Compton wavelength of the 
bosons. The equation of state ([55)1 can be written 



P = Kp^ 



with 



K 



(77) 



This returns the equation of state obtained by Arbey 
et al. [1^ for a self-interacting scalar field in the non- 



relativistic regime showing that the relation between a 
and A given by Eq. ([76]) is correct. We note that the 
parameter A can be expressed as 



A ^ 9.523^^. 
Stt 1 fm 2mn 



(78) 



We can then use (A, m) instead of (a, m) as independent 
physical variables. Finally, the dimensionless parameter 
K can be written 



K = 0.324 



Stt 



(79) 



We can now express the results in terms of A. The 
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scaling of the maximum mass is given by 



an equation of state 



A 1 



G 



(80) 

where Mp ^ {hc/Gy/'^ is the Planck mass. This is the 
scaling of the maximum mass obtained by Colpi et al. 

for a self-interacting scalar field. For A ~ 1, it is of 
the order of the Chandrasekhar mass ~ Mp/m? . On the 
other hand, the scaling of the minimum radius is given 

by 



\Gm? ) 



SnGc 



1/2 



1 




This is the scaling of the minimum radius given by Arbey 
et al. P5 for a self-interacting scalar field. Finally, the 
scaling of the maximum density is given by 



m c 



4^3 



Am c' 



(82) 



Now, the maximum mass obtained by Colpi et al. [2^ in 
the fully relativistic regime is 

M„,ax = 0.22 \/ = O.22V2A/, = 0.31Af„ (83) 
V 47r 

which is smaller than our previous estimate il/,„ax ~ 
0.5001M* based on a non-relativistic equation of state. 
Therefore, relativistic effects tend to reduce the maxi- 
mum mass. 

Colpi et al. [l^] showed that, in the Thomas-Fermi 
limit, the scalar field becomes equivalent to a fluid with 



P = 



36K 



12K 



1/2 



- 1 



(84) 



where K is given by Eq. ([77|. For p — (low or mod- 
erate densities), we recover the polytropic equation of 
state p = Kp^ corresponding to a non-relativistic BEC 
with short-range interactions. For p — > +00 (extremely 
high densities), we obtain the ultra- relativistic equation 
of state p ~ p(? jZ, similar to the one describing the core 
of neutron stars modeled by the ideal Fermi gas [13- 
[39I . We know that a linear equation of state p = qpc^ 
yields damped oscillations of the mass-central density re- 
latioii, and a spiral structure of the mass-radius rela- 
tion [39|, similarly to the isothermal equation of state 
in Newtonian gravity [13] . Therefore, our BEC model 
will exhibit this behaviour, just like standard neutron 
stars. However, our BEC model differs from standard 
neutron star models in that, at low or moderate densi- 
ties, p = Kp^ with K = 2TxaT? jir? instead oip = K' p^/^ 
with K' = (l/5)(3/87r)2/3/i/m^/^). This implies, in par- 
ticular, the existence of a maximum radius given by Eq. 
([75|) . corresponding to the Newtonian limit. 

Substituting the equation of state (|84p in the TOV 
equations, using 



P'{P) = 



1 



y/1 + UKp/c 



(85) 



and introducing the same notations as before, we obtain 



de 



6[^(Vl + 12cr6l-l)2 + cr6l] W+£:(Vl + 12cr6l-l) 



^2(1 - Aav/i){l - 1/V1 + 12(t6') 



(86) 



dv 



instead of Eqs. If we expand the square roots 

for (T ^ 1, we recover Eqs. (|58p - ([59|) . However, this is not 
a uniform expansion and the two equations ()86|) -([87 |) and 
([58ll - (j59p arc in fact different even for small values of a 
(of course, they both reduce to the Lane-Emden equation 
(IMl) for cr = 0). 

The velocity of sound at the center of the configuration 

is 



1 - 



(87) 



I 

and we always have (Cs)o < c. The series of equilibria 
becomes unstable after the first mass peak. We find that 
instability occurs for a > <t'^ with 

cr^ = 0.398. (89) 

The values of ^1 and w(^i) at this point are 

CJ = 1.923, w(^;)= 0.3865 (90) 

The corresponding values of the radius, mass and central 
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density are 



1.923 



Om- 



1/2 



4.047 k km, (91) 



M' - 307 ^ 



0.436 kMq, 



(92) 



and 



= 1.929 X 10^^ g/cm-^, (93) 



respectively. The maximum mass A/max = 0.307Af* is 
very close to the one [see Eq. (|83| ] found by Colpi et 
at. [m by solving the Klein-Gordon-Einstcin equations. 
This shows the accuracy of the hydrodynamical approach 
in the TF limit. The dimensionlcss curves giving the 
mass-central density, radius-central density, mass-radius 
relations and some density profiles are plotted in Figs.[7]- 

m 




1000 10000 



FIG. 7: Dimensionless mass-central density relation of a 
relativistic BEG with short-range interactions modeled by 
the equation of state (|84p . There exists a maximum mass 
M^^^/Mt, = 0.3076 at which the series of equilibria becomes 
dynamically unstable. The velocity of sound is always smaller 
than the velocity of light. We note that the mass-central den- 
sity relation presents damped oscillations at high densities 
similarly to standard neutron stars [3^. [3^. 



Remark: At T 
takes the form 



0, the first law of thermodynamics 



dp - 



-dn. 



(94) 



where nm is the rest-mass density. Integrating this re- 
lation with the equation of state ([84)) . we can obtain 
the relation n{p). For p — > (non-relativistic limit), 
we get p = nm, leading to P ~ Kp^ = K{nm)^ 
corresponding to a polytrope n = 1. For p +oo 
(ultra- relativistic limit), we get p oc (nm)^^'^, leading to 
P ^ pc^/3 cx (nm)^/"^ corresponding to a polytrope n = 3 




FIG. 8: Dimensionless radius-central density relation of a rel- 
ativistic BEG with short-range interactions modeled by the 
equation of state (|84p . 
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FIG. 9: Dimensionless mass-radius relation of a relativistic 
BEG with short-range interactions modeled by the equation 
of state (|84p . The mass-radius relation presents a snail-like 
structure (spiral) at high densities similarly to standard neu- 
tron stars [33,|4j|. There exists a maximum mass Mmax/M* = 
0.307 and a minimum radius Rmm/ R* ~ 1.923 corresponding 
to a maximum central density (/9c)max = 0.398p*. There also 
exists a maximum radius Rmmx/R* = tt corresponding to the 
Newtonian limit a — J> 0. 



like for an ultra-relativistic Fermi gas at T = (standard 
neutron star) . These results are consistent with those ob- 
tained by Goodman . The proper number of particles 
is 



N 



n(r) 



2GM{r) 



-1-1/2 



47rr dr. 



(95) 



It can be shown that a general relativistic, spherically 
symmetric, gaseous star at T = is dynamically stable 
with respect to the Einstein equations if, and only if, it 
is a maximum of N[p\ at fixed mass M[p\ = M (see, e.g., 
[i^li^l). The first order variations 5N—a5M = 0, where 
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FIG. 10: Dimensionless density profiles corresponding to a 
(Newtonian) and a = a'^ = 0.398 (maximum mass). 



a is a Lagrange multiplier, yield the TOV equations de- 
termining steady state solutions. The ensemble of these 
solutions forms the series of equilibria. Then, using the 
Poineare theorem [s^, one can conclude that the series 
of equilibria becomes unstable at the first mass peak and 
that a new mode of instability appears at each turning 
point of mass in the series of equilibria (see [l^l for an 
alternative derivation of these results). At these points, 
we have (5iV = bM = so that the curve N{M) presents 
cusps at each point where M{pc) reaches an extremal 
value. An illustration of this behavior is given in Figure 
5 of [1^ . These results of dynamical stability for general 
relativistic stars are similar to results of dynamical and 
thermod yna mical stability for Newtonian self-gravitating 
systems jii, 113, Ei. 



VII. ASTROPHYSICAL IMPLICATIONS 

One of the most important results in general rela- 
tivistic astrophysics is the existence of a maximum mass 
of the neutron stars [s^l- Ultra-dense compact objects 
may have a stable equilibrium configuration until their 
mass M is equal to the maximum mass Af^ax- By 
integrating the mass continuity equation and the hy- 
drostatic equilibrium equation for a star made of free, 
non-interacting, neutron gas, Oppenheimer and Volkoff 
[37I have shown that the maximum equilibrium mass is 
Mov = O.7M0, with a corresponding radius of the or- 
der of Row = 9.6 km, and a central density of the order 
of pc = 5 X lO^^g/cm'^. Using a variational method in 
which the equation of state was constrained to have sub- 
luminal sound velocity and to be stable against micro- 
scopic collapse, Rhoades and Ruffini [i^, proved that, in 
the regions where it is uncertain, the equation of state 
that produces the maximum neutron star mass is the 
one for which the sound speed is equal to the speed of 
light, i.e. p = p(? . As a result, they found a maxi- 



mum neutron star mass Mmax ~ 3.2M0, assuming un- 
certainty in the equation of state above a fiducial density 
po ~ 4.6 X 10^^ g/cm^. More realistic models that take 
into account the composition of the star and the inter- 
action between neutrons led to values of the maximum 
mass of the neutron stars in the range 1.5 — Z-IMq j47l |. 
The main reason for the lack of a better theoretical value 
of the maximum mass of the neutron star is the poor 
knowledge of the equation of state of hadronic matter at 
high densities. 

With the use of Eqs. ([Mll-dMl) it follows that the scaled 
mass, radius and central density satisfy the relations 



— = — , p^Mi 



27rG3' 



P*Rt 



2-kG 



(96) 



With the use of Eqs. dH]), dTO]), ([TTI) and (|96 l) . we ob- 
tain the following radius-mass, central density-mass, and 
central density-radius relations for the maximally stable 
non-relativistic BEG configuration. 



GM„ 



2aew(6) 



= 5.599- 



Mr. 



km, (97) 



4fT>2 (Ci) 



2^G3Af2ax 



1.026 X 10 



16 ( ^^O 



cm'^ 



r, (98) 



I-kGRI.;^ 



= 3.215 X 10^^ fl^) ^. (99) 

V -Rmin / cm-' 

For the relativistic BEG star, using Eqs. ([M]) and ([gO]) . 
we find 



i?,ni„ = 9.271%^ km, 



(100) 



(..)_. 3.682 X 10- (^) ^, (101) 



3.160 X 10^ 



15 / 10 km\ g 



cm' 



3 ■ 



(102) 



In the non-relativistic case, the mass-radius ratio of 
the star can be expressed as 



2C?A/n,ax ^ 4g,i>(6) 



0.529, 



while, for the relativistic case, we obtain 

=0.319. 



(103) 



(104) 
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A classical result by Buchdahl [48| shows that for static 
solutions of the spherically symmetric Einstein-matter 
systems, the total mass M and the area radius R of the 
boundary of the body obey the relation 2GM/c^R < 
8/9 = 0.888, the equality sign corresponding to constant 
density stars. For BEC stars, Eqs. (fTOS]) and (fTOi]) obvi- 
ously satisfy the Buchdahl inequality for the mass-radius 
ratio. 

We emphasize that the radius i?niin and the central 
density (pc)max given by Eqs. (|M)) - (jl02p only depend on 
the mass Mmax- In particular, they do not explicitly de- 
pend on the two physical parameters of the model, the 
scattering length a, and the particle mass m. On the 
other hand, the maximum mass Mmax depends on these 
two parameters only through their ratio a/m^, or equiv- 
alently, through the parameter k, and it can be obtained 
from the relations 



Mmax =O.71«;M0, 



M' 

max 



0.4368 kA/q, (105) 



respectively. Hence, all physical parameters of the model 
are determined by the mass Mmax of the star, which 
can be obtained from observations. Therefore, in the 
present model we have only one free parameter n = 
(a/fm)^^'^ (2m„/m)'^/^. With respect to a scaling of the 
scattering length and of the particle mass of the form 

a/fm — > /?! (a/fm) , TO/2m„ — > /32 (m/2m„) , (106) 

where /3i,/32 are constants, the parameter k scales as 



^ ol/2o-3/2 



bmce K = 0.324 (A/87r)^/^ (2TO„/m)2, we equivalently 
conclude that the maximum mass Afniax depends on the 
two parameters (A,m) only through their ratio A/m^. 

By assuming that the mass of the star is Af^ax = ^Mq , 
in the non-relativistic case we obtain for the param- 
eters of the star (/Oc)max = 0.256 x 10^^ g/cm'^ and 
-Rmin = 11-2 km, respectively, independently on the 
values of a and m. On the other hand in this case 
K ~ 2.816. If we take m — 2m„, this corresponds to 
a scattering length a = 7.93 fm, and a coupling constant 
A ~ 1.90 X 10"^. For the relativistic equation of state, 
we find (/9c)max = 0.091 X 10^^ g/cm^ and i?,„in = 18.54 
km. For the parameter n we obtain k = 4.578. If we 
take m = 2m„, this corresponds to a scattering length 
a = 21.0 fm and a coupling constant A = 5.02 x 10'^. 
Therefore if wc assume values of k, of the order of /c ^ 3 
in the non-relativistic regime, and k ~ 5 in the rela- 
tivistic regime, we obtain stellar objects with physical 
parameters in the range M ~ 2AfQ, i? ~ 10 — 20 km, 
and pc ~ 0.3 — 0.1 x 10^^ g/cm'^, respectively. The only 
free parameter in the model, k, uniquely determines the 
mass Afmax of the star (or conversely) . 

It may be of interest to make a connection with the re- 
sults of Oppenheimer & Volkoff [s^l- In their model, the 
mass, the radius, and the central density of the critical 



(107) 



configuration are 



A/ov = 0.376 



(Pc) 



ov 



Roy = 9.36 



3.92 X 10" 



3/2 



GMov 



= 0.7 Mr: 



= 9.6 km, 



(108) 



(109) 



5 X 10^'^g/cm^ (110) 



In our model, introducing the parameter A, the maximum 
mass is given by 



. . X f he 
A/max = 2acV(^i)\l -^{-Q 



If we write the boson mass as m 
maximum mass in the form 



3/2 



1 



(111) 



fcm„, we obtain the 



Mr, 



(112) 



where C = 0.265 for a non-relativistic equation of state 
and G = 0.163 for a relativistic equation of state. Us- 
ing the previous relations, we can then easily relate i?inin 
and (/9c)max to Roy and (pc)ov- Equation (|112p clearly 
shows that, with respect to the standard Oppenhcimer- 
Volkoff model, we have an additional parameter A (the 
strength of the self-interaction), which gives the possi- 
bility of obtaining higher values for the maximum mass. 
Using Eq. ([79l) which becomes k = we can 

rewrite Eq. (|112p as 



M^niax = 3.87 CkA/, 



ov- 



(113) 



Presently, there is conclusive observational evidence 
from pulsar studies for the existence of neutron stars 
with masses significantly greater than I.BMq [i^. By 
using the Shapiro time delay to measure the inclination, 
the mass of PSR J1614-223048 was recently determined 
to be 1.97 ± O.O4Af0 Moreover, a number of X-ray 
binaries seem to contain high-mass neutron stars: about 
1.9Af0 in the case of Vela X-1 and 2AMq in the case 
of 4U 1700-377 [4^. Even more intriguing is the case of 
the black widow pulsar B1957-I-20, with a best mass esti- 
mate of about 2.4Af0 [5l[. This system has both pulsar 
timing and optical light curve information. B1957-I-20 is 
located in an eclipsing binary system, consisting of the 
1.6 ms pulsar in a nearly circular 9.17 h period orbit, and 
an extremely low mass companion, A/c w O.O3A/0. It is 
believed that irradiation of the companion by the pulsar 
strongly heats its cosmic environment to the point of ab- 
lation, leading to a comet-like tail, and a large cloud of 
plasma. The plasma cloud is responsible for the eclips- 
ing. The pulsar is literally consuming its companion, 
hence the name black widow. The mass of the compan- 
ion star has been reduced to a small fraction of its original 
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mass [4^. On the other hand, a measured mass of 2AMq 
would be incompatible with hybrid star models contain- 
ing significant proportions of exotic matter in the form 
of hyperons, some forms of Bose condensates, or quark 
matter [H]. 

However, the mass and radius of the 2 — 2AMq neu- 
tron stars perfectly fit the expected properties of a Bose- 
Einstein condensate star. For k ~ 3, the mass of a typ- 
ical general relativistic Bose-Einstein condensate star is 
of the order of two solar masses, with a radius of around 
11 km. Therefore, we propose that the recently observed 
2 — 2AMq mass neutron stars could be typical Bose- 
Einstein condensate stars. 

A last comment may be in order. If we apply the 
same model (self-gravitating BEG with short-range in- 
teractions) to dark matter |26| . and use the Newtonian 
approximation (which is valid in this context), the radius 
of a dark matter halo is given by Eq. (|33|) . which can be 
rewritten 

i? = 1.746 X 10- (^) kpc. 

(114) 

Again, we note that the radius R determines the ra- 
tio ajrr? or Xj'nv' . Estimating the radius of dark mat- 
ter halos by i? = 10 kpc, we obtain vn? ja = 3.049 x 
lO-*' (eV/c2)3/fm and mVA = 23.94 (eV/c^)* 

VIII. DISCUSSIONS AND FINAL REMARKS 

In the present paper, we have proposed that the core of 
neutron stars is a superfluid in which the neutrons form 
equivalent of Cooper pairs, so that they act as bosons of 
mass 2m„. Therefore, once the Bose-Einstein condensa- 
tion takes place, the neutron star should be modeled as a 
self-gravitating BEG star (boson star). In our approach, 
we also assume that the bosons have a self-interaction, 
described by a scattering length a. The basic properties 
of the gravitationally bounded Bose-Einstein condensates 
have been obtained in both Newtonian and general rela- 
tivistic regimes. To obtain the physical characteristics of 
the system, we have used the Madelung representation 
in which condensates can be modeled by using the hy- 
drodynamic Euler equations describing a gas whose den- 
sity and pressure are related by a barotropic equation of 
state. For the study of the Bose-Einstein condensate we 
have adopted the Thomas-Fermi approximation, which 
is valid if the total number of particles N obeys the con- 
dition N ^ R/ira [l^, H^, Ho], a condition which can 
be reformulated, with the use of the mass density of the 
condensate, as p ^ 'im/ Air^R^a. This restriction is obvi- 
ously satisfied by condensates with densities of the same 
order as the nuclear density. 

In the physically most interesting case, corresponding 
to a quartic non-linearity term in the energy functional, 
the equation of state of the Bose-Einstein condensate is 



that of a polytrope with polytropic index n = 1. In 
this case, the radius and the mass of the Newtonian stel- 
lar condensate can be obtained in an exact form. In 
a Bose-Einstein condensed neutron star the mass m of 
the particle does not need to coincide with the neutron 
mass. For the mass of the condensed particle we have 
used an effective value of the order of m* = 2to„. This 
value is justified by the high densities in the neutron star 
cores, where the process of Bose-Einstein condensation is 
most likely to occur via the formation of Cooper pairs. 
However, we have also explicitly presented the numerical 
values of the basic physical parameters of the stars for 
other values of the mass. 

General relativistic effects impose strong constraints 
on the maximum mass. In the framework of the general 
relativistic approach one must numerically integrate the 
structure equations of the star. In this way, we obtain 
a large class of stable astrophysical objects, whose ba- 
sic parameters (mass and radius) depend on the particle 
mass m and scattering length a. Since the values of a and 
m are not well-known, this offers the possibility to obtain 
a maximum mass for neutron stars that is larger than the 
Oppenheimer-Volkoff limit of 0.7 M©, and may be com- 
patible with recent observational determinations of the 
masses of some neutron stars. This is possible because 
we have two new parameters in our model, the boson 
mass m and the scattering length a, which give additional 
freedom (although these parameters should be ultimately 
determined by fundamental physics). We have found 
that the maximum mass Afniax('*) of the condensate star, 
given by Eqs. ([75]) or ([M)) . depends in fact on a single pa- 
rameter K{a,m) which is proportional to the ratio a/m^ 
(or, equivalently, to the ratio \/m^). Since the radius 
Rmin and the density (pc)max only depend on Mmax, all 
the physical properties of the BEG stars are determined 
by the parameter k. Condensates with particle masses of 
the order of two neutron masses and scattering length of 
the order of 10 — 20 fm (corresponding to k ^ 3 — 5) have 
maximum masses of the order of 2Mq, minimum radii 
in the range of 10 — 20 km and maximum central den- 
sity pc ^ 0.3 — 0.1 lO^^g/cm'^ in the non-rclativistic and 
relativistic regimes, respectively. On the other hand, for 
a = 1 fm, the maximum mass of the condensate varies 
between 0.4 — 0.7 for m = 2to„ (corresponding to 
K = 1), between 1 — 2 Mq for m ~ m,n (corresponding to 
K = 2.8), and between 10 — 16 M0 for m = to„/4 (corre- 
sponding to K = 22), a value which may correspond, for 
example, to the effective (density dependent) kaon mass 
to|^ in the interior of neutron stars. Kaon condensation 
may provide an important example of Bose-Einstein type 
stellar condensate [l3|. Attraction from nuclear matter 
could bring down the mass of the kaon to an effective 
value of mj^, « 200MeV/c2 « to„/9.38 MeV/c^. From 
Eq. ([55]) it follows that Newtonian kaon condensates with 
kaon effective mass of the order of mj^ = m„/10 could 
have masses as high as lO^M©, and radii of the order of 
600 km [see Eq. ([55]) ]. Thus, Bose-Einstein condensate 
stars, formed from small mass particles, may represent vi- 
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able candidates for the super-massive "black holes" that 
reside at the galactic centers. However, general relativis- 
tic effects strongly restrict the value of the maximum 
mass of super-massive BEC stars. From Eq. ([^^ it fol- 
lows that for m = mj^ = m„/10 we have n = 89.44. 
With the use of Eq. (|73p we obtain a maximum mass 
of M = 63.50 A/q for the kaon condensate star, with a 
radius R = 355 km [see Eq. ([7^ ]. On the other hand, 
smaller mass condensed particles can have significantly 
higher maximum relativistic masses. Hence the Bose- 
Einstein condensation process in the early universe may 
have provided the seeds from which super-massive black 
holes were eventually formed through accretion of inter- 
stellar matter. 

Presently, the mass of the neutron stars can be deter- 
mined very accurately, and many of them have masses 
in the range of 2 — 2.4 solar masses, which arc very diffi- 
cult to explain by the standard neutron matter models, 
including those with exotic matter like quarks. However, 
these mass values could be very easily explained by our 
model if neutron stars can be modeled as BEC stars. 

Bose-Einstein condensate stars could have a normal 
matter crust, since we expect that the condensation can- 
not take place at densities smaller than the nuclear den- 
sity or quark deconfinement density. The presence of 
the thin crust increases the mass and the radius of the 
condensate star by a factor of 10% or 17%, respectively. 
Therefore, the presence of a neutron crust does not mod- 
ify significantly the basic physical properties of the star. 
Distinguishing between Bose-Einstein condensate stars 
and "standard" neutron stars or other type of conden- 
sate or quark stars could be an extremely difficult obser- 
vational task. Similarly to the case of quark stars [s^lssj . 
we suggest that high energy radiation processes from the 
surface of the condensate may provide some distinctive 
features allowing a clear differentiation of these different 
types of stellar objects. 



In a very general approach, one may assume that the 
masses m of the particles forming the stellar type con- 
densate are anisotropic, and they should be described by 
a mass tensor rriij. Such anisotropic masses are known 
from condensed matter physics where they are encoun- 
tered in effective mass calculations for electrons immersed 
in a band structure, in the case of excitons (electron-hole 
couples held together by the Coulomb attraction) and 
in BEC for semiconductors. The doping structure of the 
semiconductor and its anisotropies would give place to an 
effective mass matrix for the paraexcitons (singlet exci- 
tons) at least in the low momentum approximation Q . A 
different value for the effective mass m may considerably 
increase (or decrease) the total mass of the condensate. 

A rotating Bose-Einstein condensate may exhibit a 
very complex internal structure and dynamics, mainly 
due to the presence of vortex lattices. The vortex lat- 
tices may evolve kinetically, with each vortex following 
the streamline of a quadrupolar flow. The quadrupolar 
distortions can lead to a disordering of the vortex lat- 
tice, and to an instability due to inter-particle collisions, 
flnite temperature effects or to the quadrupolar distor- 
tions induced by the external potential. On the other 
hand, due to the high neutrino emissivity, which is sig- 
nificantly enhanced due to the condensation, kaon con- 
densate stars are very dark objects. Hence their obser- 
vational detection may prove to be an extremely difficult 
task. The possible astrophysical/observational relevance 
of these processes will be considered in a future publica- 
tion. 
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